Emergence of stability in a stochastically driven pendulum: beyond the Kapitsa effect 
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We consider a prototypical nonlinear system which can be stabilized by multiplicative noise: an 
underdamped non- linear pendulum with a stochastically vibrating pivot. A numerical solution of 
the pertinent Fokker-Planck equation shows that the upper equilibrium point of the pendulum can 
become stable even when the noise is white, and the "Kapitsa pendulum" effect is not at work. The 
stabilization occurs in a strong-noise regime where WKB approximation does not hold. 
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It has been known for a long time that multiplica- 
tive noise can enhance stability of nonlinear systems. 
Examples are numerous indeed and culminate at noise- 
induced phase transitions far from equilibrium [1]. This 
paper deals with a noise-induced stabilization of oscillat- 
ing systems. As a prototypical example we consider an 
underdampled nonlinear pendulum with a stochastically 
vibrating pivot. The stochastic driving introduces both 
multiplicative and additive noise, see Fig. [TJ Our numeri- 
cal simulations clearly show that the multiplicative noise 
can stabilize the otherwise unstable upper equilibrium 
point of the pendulum. The mechanism for this stabi- 
lization is markedly different from, and more subtle than, 
the "Kapitsa pendulum" mechanism. The Kapitsa pen- 
dulum involves a (deterministic) monochromatic para- 
metric driving of the pendulum at a frequency that is 
much higher than the natural frequency of the pendulum 
Here the upper equilibrium point becomes stable if 
the driving acceleration is higher than a critical value 
depending on the pendulum length and the gravity ac- 
celeration. In the Kapitsa pendulum problem the change 
of stability of the upper equilibrium point comes from a 
change in the effective potential of the pendulum [2]. 
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FIG. 1: (color online). Schematic of the stochastically driven 
simple gravity pendulum. 

Extensions of the Kapitsa pendulum effect to mul- 
tiplicative stochastic driving have also been considered 
[3, |4), see Ref. [5] for a review. In these extensions the 
noise spectrum is strongly peaked at a single frequency 
which is much higher than the natural frequency of the 
pendulum. The presence of a high-frequency noise of a 



sufficient strength modifies the effective potential which 
can stabilize the upper equilibrium point. Theory- wise, 
this setting introduces a time-scale separation which per- 
mits a perturbative treatment [Hj]. In this work we con- 
sider a model stochastic driving with a fiat spectrum: a 
white noise. Here all frequencies from to oo are equally 
present, and there is no time-scale separation. We will 
show that, for such a noise, the upper position of the pen- 
dulum can also become stable. However, the stabilization 
cannot be traced to a change in the effective potential of 
the pendulum. 

A stochastically driven simple gravity pendulum can 
be described by a Langevin equation: 

= fi, (1) 
n = -u)l sin - 2jQ + - sin y/JI & (t) + y/a 6 (t) , (2) 

where is the deviation angle of the pendulum, see Fig. 
[U Q is the angular velocity, uo = \fgjl is the harmonic 
frequency of the pendulum, 7 is the damping factor, I is 
the the pendulum length, g is the gravity acceleration, 
and \i and a are the magnitudes of the multiplicative 
and additive noises £1 and £2, respectively. The noises 
are assumed to be guassian, white with zero mean and 
mutually uncorrelated: 

&(*)> = 0, <&(W)> = m-t')k) , hj = 1,2. (3) 

The Langevin equations (pQ) and ([2]) are equivalent (see, 
e.g. Ref. [6[) to the following Fokker-Planck equation for 
the probability distribution W(0,£l,i): 



W t = -nWe+Uo s'mOWn 



(4) 



where the indices Q and t denote the corresponding 
partial derivatives of W{6, Cl,t). As the noises are Q- 
independent, there is no difference between the Ito and 
Stratonovich interpretations. Introducing the dimension- 
less variables t = ou^t, Vt = Q/uo and W(t, 6, Q) = 
(Jo W(t, 0, Q), we can rewrite the Fokker-Planck equation 
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in a dimensionless form: 

W t = -QWe + anOWa 
+ (s + 5sm 2 e) Wi 



(5) 



where T = 7/cJo, £ — ol/uq anc ^ ^ = M/G^o) are ^ ne 
rescaled parameters of the system, and the tildes are 
omitted. 

We assume that, after a transient, the stochastic 
system approaches a smooth steady state for which 
W(9,Q,t) is independent of time: W(6,Vt,t — > 00) = 
VK(#,^). The steady-state probability distribution 
W(9, f2) obeys the equation 



d 



-nWo + 8m6Wn 
+ (£ + (5sin 2 ^) Wh n = 0. 



(6) 



We classify a point (0, f2) as a stable point of the system 
if it is a local maximum of the stationary probability 
distribution W(9,Q). The necessary and sufficient con- 
ditions for a function of two variables f(x,y) to have a 
local maximum at (#0,2/0) are (see, e.g. Ref. Q): 

fx(x , yo) = f y (x , yo) = , (7) 
fxx(xo, yo) < or f yy (x , yo) < , (8) 

fxxfro, yo) fyy(xo, 2/0 ) /^fao, 2/0 ) > , (9) 

where the indices x and 2/ denote partial derivatives. Let 
us examine the stability properties of the upper equi- 
librium point (9 = 7r, Q = 0) of the driven pendu- 
lum. Equation (|6|) is invariant under the transformation 
— > 2tt — 0,Q — > — fi, that is under reflection of the axes 
and f} around the point (7r, 0). Its solution W"(0, f2) must 
obey the same symmetry. Therefore, the first derivatives 
We and Wq must vanish at (7r,0), and so Eqs. (0 are 
satisfied there. This immediately follows 



Whn(7r,0) 



2r 



-W(tt,0) < 0, 



so Eq. (j8j) is also satisfied at (tt,0). As a result, the nec- 
essary and sufficient condition for (tt, 0) to be a stable 
point is given by Eq. ([9]): 

A = W^n(7r,0)W-^(7r,0) - W^(tt,0) > 0. (10) 

For S = (only additive noise), the steady-state equation 
(|6]) is soluble analytically, see Ref. [8[: 



W(0,ft) = 



pl/2 



2^/% 1 /2/ ( 2 r/£) 



exp 



--to 2 



2cos6») 



(11) 

where io(- • • ) is the modified Bessel function. In this case 
the point (-zr, 0) is unstable, as the stability parameter A, 
defined in Eq. (jTUJ) , is negative: 



7r 3 e 3 I^(2T/e) 



exp 



4r A r, 

- <0. 



(12) 




FIG. 2: (color online). A steady-state probability distribution 
with no multiplicative noise. The parameters are e — 0.1, 
r = 0.1 and 5 — 0. This solution evolved from the initial 
condition W = (2ny/ : K)~ 1 exp(— Q 2 ) after 60 normalized time 
units. It coincides with the analytical solution (|11|) . 



The probability distribution ([TT]) is depicted in Fig. [2l 

What happens when S > 0? First, we show that there 
is no stabilization in the weak- noise limit, {~e<l,r. 
In this limit one can use a dissipative variant of WKB 
approximation, see e.g. Ref. @, and make the ansatz 



S(0,fi) 



W(0, SI) = A{9, n) exp 



assuming that the pre-factor A varies on scales much 
larger than \je. In the leading order in 1/e one obtains 



nS e - sin 0Sq - 2TttSn + ^1 + - sin^ Oj = . (13) 

This first-order PDE for the action S(9, Q) has the form 
of a stationary Hamilton- Jacobi equation (see e.g. Ref. 
Q) with the Hamiltonian H(9,Q,pi,p2) of the form 



H = n Pl - sm0p 2 - 2rn P2 + ^1 + - sin^ 6 J pi . (14) 

Here p\ = Sq and P2 = Set are the canonical momenta 
conjugated to the coordinates and ^, respectively. As 
the Hamilton-Jacobi equation ([T3j) is stationary, we are 
interested in the zero-energy dynamics H = 0. The 
Hamilton's equations 



= fi, ft = -sin0-2rft + 2 ^l + -sin 2 0j p 2 , (15) 

pi = cos0p2 sin 20 pi ? P2 = — Pi + 2Fp 2 , (16) 

have two zero-energy fixed points: a\ = (0, 0, 0, 0) and 
&2 = (tt, 0,0,0), corresponding to the lower and upper 
equilibrium points of the pendulum, respectively. The 
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linear stability of the upper equilibrium point a 2 is de- 
termined by the quadratic approximation to the Hamil- 
tonian around a 2l 

H(9, n,pi,p2) - Q Pl + 0p 2 - 2YQp 2 + p\ . (17) 

This quadratic approximation is noiseless, as the noise 
term in Eq. ([T4]) . (S/e) sin 2 Op?;, is bi-quadratic in small 
deviations from the fixed point. Therefore, a weak mul- 
tiplicative noise does not generate any correction to the 
potential of the pendulum, in agreement with an early 
observation [3[, and cannot change the stability proper- 
ties of the system. 

However, for a strong multiplicative noise numerical 
solutions of the time-dependent Fokker-Planck equation 
© do show emergence of stability of the upper equilib- 
rium point. We obtained these numerical results with a 
"Mathematica" PDE solver. The numerical domain was 
< 6 < 2tt and — Vtmax < ^ < Qmax with flmax chosen, 
separately for each set of parameters, sufficiently large. 
Periodic boundary conditions in were used. As to the 
boundary conditions in ^, we checked that, at sufficiently 
large Qmax, the steady-state solution remained the same 
up to 1 per cent, whether we imposed periodic or zero- 
W conditions at the ^-boundary. Larger ^t ma x needed 
to be taken when T became smaller than e and when 5 
became large compared to T and e. The values of ^l ma x 
that we used ranged from fl m ax = 2 for T = 0.5, e = 0.05 
and S ~ 0.15 to n max 15 for T 0.025, e = 0.5 and 
S ~ 0.43. After verifying that Vtmax is sufficiently large, 
we used the zero-W boundary conditions in as this 
choice reduced the computational time. For each set of 
parameters 5, Y and e we used the analytical solution ([TT]) 
for 8 = as the initial condition. 

We ran the solution until t = t max such that the value 
of A, evaluated at t = t max , was within 1 per cent from 
its value at t = t max /2. The larger the parameter T, the 
smaller t max was needed. The values of t max that we 
used ranged from 40 for Y = 1, e = 0.3 and 5 = 0.6 to 
140 for T = 0.025, e = 0.05 and S = 0.15. 

At small Y/e or large 5 the steady-state probability 
distribution broadens which demands a larger D maa: , a 
longer computation time and more computer memory. 
At large Y/e or small 5 the distribution becomes too nar- 
row to numerically resolve it with confidence. We ver- 
ified that, starting from an arbitrary initial condition, 
the probability function always converges to the same 
steady-state solution. We also checked that for S = 
the steady-state solution coincides, with a high accuracy, 
with the analytical solution ([TT]) . 

Figure [3] gives a typical numerical example of a steady- 
state probability distribution having a distinct local max- 
imum at (7r, 0), in addition to the expected (higher) peak 
at (0,0). Also noticeable is a significant broadening of 
the probability distribution, compared with the case of 
only additive noise, see Fig. [2] 




FIG. 3: (color online). Multiplicative noise stabilizes the up- 
per equilibrium point of the pendulum. Shown is the steady- 
state probability distribution W{6, fi) which exhibits a local 
maximum at the upper equilibrium point (6 = 7r, Q = 0). The 
parameters are e = 0.1, Y = 0.1 and 6 = 2. The steady-state 
probability distribution was found by solving numerically the 
time-dependent Fokker-Planck equation J5]). It evolved from 
the initial condition (|11[) after 60 dimensionless time units. 

To determine the stability properties of the upper equi- 
librium point (7r,0), we plotted the stability parameter 
A, defined in Eq. ([TO]), versus the rescaled magnitude of 
the multiplicative noise J, for different values of Y and e. 
Several examples of such plots are shown in Fig.[H Using 
interpolation, we found the critical values 5 C = 5 c (Y,e) 
of (5, at which A = 0. 
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FIG. 4: (color online). The stability parameter A, normalized 
to |Ao| from Eq. (|12|) . is plotted versus the rescaled magnitude 
of the multiplicative noise 6 for different values of e at Y = 0.1. 
The lines are only given for guiding the eye. The crossing 
point of each plot with the line A = yields S c . 

Figure [5] shows the plot of 5 C versus Y at different e. 
One can see that S c slowly decreases with Y and ap- 
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FIG. 5: (color online). The critical rescaled magnitude of the 
multiplicative noise 5 C vs. the rescaled damping factor Y for 
different values of e. The lines are only given for guiding the 
eye. 
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FIG. 6: (color online). The critical rescaled magnitude of 
the multiplicative noise S c vs. the rescaled magnitude of the 
additive noise e for different values of Y. The lines are only 
given for guiding the eye. 



proaches a finite value at Y — > 0. The plot of S c versus e at 
different T, presented in in Fig.[6l shows a much stronger, 
square-root-like dependence at small e. Our data for dif- 
ferent T do not contradict a power-law behavior S c ~ e a 
with a ~ 0.44, as shown in Fig. [Tlfor Y = 0.1. When e 
approaches 0, S c also goes to zero, and one would expect 
to always see stability in this case. As e — > 0, however, 
the probability distribution develops singularities both 
at (0,0), and at (ir,0), but the peak at (0,0) becomes 
much higher than that at (7r, 0). A local maximum at 
(7r, 0) is physically insignificant in this case. Only with 



a sufficiently large additive noise the probability distri- 
bution becomes sufficiently broad to allow a reasonable 
probability for the pendulum to be at and around (7r,0). 

In summary, we have investigated the stabilization of a 




FIG. 7: The log-log plot of S c vs. e for Y = 0.1. A linear fit 
gives a slope of 0.44±0.01 with the coefficient of determination 
R 2 = 0.99726. 



prototypical nonlinear oscillating system by a multiplica- 
tive white noise. The stabilization is clearly observed 
in the numerical solution of the Fokker-Planck equation 
and requires a super-critical noise magnitude. The sta- 
bilization cannot be traced to a change in the effective 
potential of the system and is not predicted by a WKB 
analysis which assumes a weak noise. 
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